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We consider a Kerr black hole acting as a gravitational deflector within the geometrical optics, 
and point source approximations. The Kerr black hole gravitational lens geometry consisting of an 
observer and a source located far away and placed at arbitrary inclinations with respect to the black 
hole's equatorial plane is studied in the strong field regime. For this geometry the null geodesies 
equations of our interest can go around the black hole several times before reaching the observer. 
Such photon trajectories are written in terms of the angular positions in the observer's sky and 
therefore become "lens equations". As a consequence, we found for any image a simple classification 
scheme based in two integers numbers: the number of turning points in the polar coordinate 0, 
and the number of windings around the black hole's rotation axis. As an application, and to make 
contact with the literature, we consider a supermassive Kerr black hole at the Galactic center as 
a gravitational defiector. In this case, we show that our proposed computational scheme works 
successfully by computing the positions and magnifications of the relativistic images for different 
source-observer geometries. In fact, it is shown that our general procedure and results for the 
positions and magnifications of the images off the black hole's equatorial plane, reduce and agree 
with well known cases found in the literature. 
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I. INTRODUCTION 

As is well known, the deflection of light by a gravitating 
body was one of the first predictions of Einstein's Gen- 
eral Theory of Relativity to be observationally confirmed. 
Later on, Einstein himself predicted what is called to- 
day a microlens: the momentarily increase in apparent 
brightness of a background star as it passes close to a fore- 
ground massive body Ij. Both the deflection of light, and 
the change in apparent brightness of a radiation source 
by an external gravitational field, are collectively known 
as a gravitational lens. 

Nowadays, gravitational lensing is a very active area of 
research, and it has found applications ranging from the 
search of extrasolar planets and compact dark matter to 
estimate the value of the cosmological parameters |2, IM 
1^. In most of these applications it is only necessary to 
assume that the gravitational field is weak and that the 
deflection angle (due to a spherically symmetric body 



of mass M) can be approximated by: 



AGM/bc 



where b is the impact parameter [^. On the other hand, 
it is well known that for a Schwarzschild black hole, the 
deflection angle diverges as 6 — > i^/iM, allowing photons 
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to orbit the black hole many times before reaching the 
observer. This gives rise to an infinite set of images at 
both sides of the black hole 0, Q- Notice that in this 
region the gravitational field is no longer weak, and the 
above approximation fails. 

Recently, a paper by Virbhadra and Ellis have renewed 
interest in such images, which they called relativistic im- 
ages 1^. Later on, Bozza 0, ^3 and Eiroa, Romero 
and Torres developed an approximation method for 
the case of strong spherically- symmetric gravitational 
fields. In fact, by expanding the deflection angle near 
the point of divergence, these authors were able to find 
analytic expressions for the positions and magnifications 
of the resulting relativistic images. Interestingly, such 
images were only characterized by the number of wind- 
ings around the black hole [23 |. 

Although several authors have began to study gravita- 
tional lenses with a rotating black hole as a gravitational 
deflector (see and references therein) , most of their 
approaches to this subject still use the weak fleld ap- 
proximation, and focus only in null geodesic motion at 
the black hole's equatorial plane. In passing, we shall 
mention for the interested reader a nice discussion of a 
Kerr black hole as a gravitational lens by Bray . Un- 
fortunately, in this work there are approximations that 
are valid only for small deviations from the straight line 
path, and therefore are not suited for studying relativis- 
tic images. If we want to consider the phenomenology of 
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the relativistic images in the strong gravitational field of 
a Kerr black hole, no approximations can be taken and 
we need to work with the full equations of motion for null 
rays. 

There have been numerous articles about the motion of 
null rays in the gravitational field of a Kerr black hole 0, 
El El 113 ■ Although some of them have address the 
gravitational lens problem, most have concentrated in 
the observational effects on accretion disks and sources 
orbiting at the equatorial plane of the black hole and 
have not discussed the phenomenology of the relativistic 
images. 

Recently, Bozza E9 studied "quasi-equatorial" orbits 
of photons around a Kerr black hole and provided an- 
alytical expressions for the positions and magnifications 
of the relativistic images. However, these approximations 
fail when the observer and the source are located far away 
from the equatorial plane. Moreover, as we show in Sec. 
V, Bozza's procedure begins to fail when the black hole's 
angular momentum increases near its maximum value, 
even if the observer and the source are close to the rotat- 
ing black hole's equatorial plane. 

In this paper, we shall discuss the phenomenology of 
the relativistic images by using the exact null equations 
of motion, with the assumption that the observer and 
the source are far away from the black hole. Although 
this new procedure is far more complicated than previous 
works on this subject, it allow us to calculate, estimate, 
and discuss for the first time the observational proper- 
ties of the relativistic images for arbitrary source and 
observer inclinations in a Kerr gravitational lens. 

Therefore the purpose of this paper is two-fold. First, 
to extend the study of relativistic images for the case 
when the gravitational defiector is a rotating Kerr black 
hole. Secondly, since the trajectory of a photon will not 
always will be confined to a plane, we shall also be con- 
cerned with photon trajectories off the black hole's equa- 
torial plane. As a consequence of this undertaking, we 
show below that all relativistic images deflected by a ro- 
tating Kerr black hole are characterized by only two in- 
tegers numbers: namely, the number of turning points 
in the polar coordinate 9, and the number of windings 
around the black hole's rotation axis. 

To facilitate reading this paper is divided as follows: in 
Sec. II to make this study self-contained a review of the 
null geodesic's equations of motion in the Kerr space-time 
is briefly discussed . In Sec. Ill we explain the gravita- 
tional lens geometry, and present a general classiflcation 
for all images that allow us to define formally what we 
mean by "relativistic images" . In Sec. IV, analytical ex- 
pressions for the images' magnification arc derived. In 
Sec. V, we calculate the positions and magnifications of 
the relativistic images for several special cases, and com- 
pare our results with those in the literature. Finally, in 
Sec. VI a discussion of our results is undertaken. An 
appendix is also included to show how the null equations 
of motion can be solved in terms of elliptic integrals. 



II. NULL GEODESICS IN THE KERR 
SPACE-TIME 



In this section we briefly review the equations of mo- 
tion for light rays in the Kerr space-time. We use the 
usual Boyer-Lidquist coordinates which, at infinity, are 
equivalent to the standard spherical coordinates. We also 
discuss about the relevant range of coordinates and con- 
stants of motion for the case of an observer and a source 
located far away from the black hole. Details about the 
relation between the null geodesic equations, the gravita- 
tional lens geometry, and observable quantities are given 
in the Sec. III. 

By a convenient choice of the afflne parameter, null 
geodesies in Kerr space-time can be described by the fol- 
lowing first-order differential system [l9|: 
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(1) 
(2) 

(3) 

(4) 



where 



R{r) = + {a^ - - 

+2 [{a -Xf+T]]r- o?ri 

6(6') = ?7 -I- cos^ 6* - cot^ 61 , 

P = (r2 + a2^_Aa, 

E = r"^ + cos^ 6 , 

A = - 2r -f . 



(5) 
(6) 
(7) 
(8) 
(9) 



In the above, = dx^/dr is the four- velocity (r is an 
affine parameter), A and rj are the constants of motion, 
a — J/M is the black hole's angular momentum per unit 
mass, and units are chosen such as GM/c^ = 1. As is 
well known, for a Kerr black hole, the parameter a is 
restricted to < |a| < 1. From Eqs. Q - igj it follows 
that the relevant integrals of motion are 
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dr 



A0 
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a{2r — aX) 



d0 



(10) 



dr 



dO 



'•r^(r^ + a^) + 2Qr(a-A) 
±A[i?(r)]' 
" a2 cos2 



(11) 



(12) 



The signs of ^/ R{r) and ^/Q{e) are those of V and C/^ 
respectively. Thus, the positive sign is chosen when the 
lower integration Hmit is smaller than the upper limit, 
and the negative sign otherwise. 

For an observer and a source located far away from the 
black hole, the relevant radial integrals can be written as 
follows: 



(13) 



we see that r] could be negative. However, assuming 
a photon crosses the equator {9 = 7r/2), implies that 
V ~ Ug > 0. Since in this research we consider a source 
behind the black hole and we are mostly interested in 
images formed by photons that go around the black hole 
before reaching the observer, we will only consider the 
case of positive rj. 

For a photon to be able to return to infinity, we need 
dW/dr > at r r,„i„. Setting W = = dW/dr we 
get: 
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r2(r-3) + a2(r + l) 
a(l -Zl) 

[4a2 - r{r - if] 



a2(l-r)2 



(17) 
(18) 



where r is the lower bound of 



with r/i < Zl < oOj 



and Th — ^ -\- vT— is the Kerr black hole's horizon. 
This is a parametric curve in the (A, rj) space and it is 
shown in Fig. ^ for the case a = 0.9. Photons with con- 
stants of motion inside the shaded region do not have a 
turning point outside the horizon so they will fall into 
the Kerr black hole. This is analogous to the motion of 



where rmin is the only turning point in the photon's tra- 
jectory, and it is defined by the largest positive root of 
R{r) — 0. For the angular integrals, however, there could 
be more than one turning point. In consequence, the 
most general trajectory is described by 



/"^min / max Z*^ 



■ ± 



(14) 



The turning points in 9 arc defined by Q{9) — and are 
given by 



cos (^min / max) 



a-' 
A2 



-I 1/2^ 
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a'' 



(15) 



It is easy to show from Eq. (jSJ that for any pair of 
parameters (A, i]), the motion in 9 is bounded by 0max and 
6'niin- Thus, the space of parameters is restricted because 
^min < (^ot9s < ^max- This constrain will be implemented 
in Sec. Ill after discussing the lens geometry. Finally, the 
deflection angle Acf) has to be chosen to satisfy a given 
source-observer geometry (see Sec. Ill for details). 

Now we want to know what region of the parameter 
space (A, 77) correspond to photons that after reaching 
''mill can escape to infinity. Writing Eq. ||2J) as 




cos^ 9 + X^ cot^ 



(16) 



FIG. 1: Photons with parameters (A,r;) inside the shaded 
region will fall into the black hole. The parametric curve is 
for a = 0.9 and is defined by Eqs. ([T71 and llTSll . 

null geodesies around a Schwarzschild black hole, where 
any photon with impact parameter h < 3\/3 can not 
escape to infinity ^ 0. This value of the parameter 
6, correspond to rmin = 3 which defines the "photon 
sphere". In the parameter space (A, 77), this correspond 
to a close region bounded by the line 77 = and the curve, 
\/A2 + T] = 3^3. This forbidden region in the parameter 
space shall be called: photon region. The radial integrals 
of Eqs. l(Tnil and ((TT|l diverge at the boundary of the pho- 
ton region (except at the line r] = Q) and take complex 
value inside of it. Therefore, knowledge of the mapping 
of this region will allow us to ensure that the equations 
of motion will remain well behaved. 
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III. THE KERR GRAVITATIONAL LENS 
GEOMETRY 



Changing to spherical coordinates and using the equa- 
tions of the straight hne, is easy to show that 



Since Kerr space-time is asymptotically flat, an ob- 
server far away from the black hole (ro ^ 1) can set 
up a reference euclidean coordinate system (x, y, z) with 
the black hole at the origin (see Fig. [SJ. The Boyer- 
Lidquist coordinates coincide with this reference frame 
only for large r. The coordinate system is chosen so that, 
as seen from infinity, the black hole is rotating around 
the z axis. For a > the rotation will be assumed to 
be in the counterclockwise direction as seen from the 
positive z axis. Without loss of generality, we choose 
(j)o — so that the coordinates of the observers in the 
Boyer-Lidquist system are {ro,0o,0). Similarly, for the 
source we have {rs,9s, (j^s)- In the observer's reference 



a, = -r„ smt 



^3^ 



dr 

2 d9 
r — 
° dr 



(19) 
(20) 



By using Eqs. ITJ - Q in lO and and further 
assuming oo, it is possible to relate the constants 

of motion A and rj to the position of the images in the 
observer's sky: 



obs&vei' 




-ai sm( 



sowre 



FIG. 2: Geometry of the rotating gravitational lens. An ob- 
server far away from the black hole, can set up a reference co- 
ordinate system {x, y, z) with the black hole at the origin. The 
Boyer-Lidquist coordinates coincide with this system only at 
infinity. The reference frame is chosen so that, as seen from 
infinity, the black hole is rotating around the z axis. In this 
system, the line joining the origin with the observer is normal 
to the Q-/3 plane. The tangent vector to an incoming light ray 
defines a straight line, which intersects the plane at the 
point (Qi,A). 



77 « (a- - a^) cos^ Oo + Pf 



(21) 
(22) 



These equations can, in turn, be written in terms of the 
angles {xi,yi) in the observer's sky by: ai w roXi, f3i w 
"ToUi- We would like to mention that Eqs. H19|l - H22|l are 
not equivalent to Eqs. (7), (8), (11) and (12) of Ref . Q^. 
However, our equations exactly coincide with Eqs . (28a) 
and (28b) of Ref. We believe that in Ref. there 
is a calculation mistake in the intersection of the tangent 
to the light ray and the a-/3 plane. The author in Ref. 
[T^ also ignore any contribution of the black hole's spin 
when expanding (d9/dr)\r^ and (d4i/dr)\r^. 

As can be seen later, it is also useful for our purposes 
to write the source's angular coordinates {0s, 4>s) in terms 
of its position in the observer's sky. From the geometry 
of Fig. 12 we obtain 



r, smt/, cos( 



f o — rs (sin 60 sin 9s cos 4>s + cos 9o cos 9s 



(23) 



^ (sin 9o cos 9s — cos 9o sin 9s cos 4>s) ^^^-j 
fo — fs (sin 9o sin 9s cos (ps + cos Bo cos 9s) ' 



frame, an incoming light ray is described by a paramet- 
ric curve x{r), y{r) , z{r), where = -f- + z^. 
For large r this is just the usual radial coordinate in the 
Boyer-Lidquist system. At the location of the observer, 
the tangent vector to the parametric curve is given by: 
{dx/dr)\ X -I- {dy/dr)\^^ y -I- {dz/dr)\ z . This vector 
describes a straight line which intersect the a-(3 plane 
shown in Fig.[5]at (a^, (3i). A line joining the origin with 
the observer is normal to this plane. We call this plane 
the observer's sky. The point (at, Pi) in this plane is 
the point {—Picos9o,ai, PiSm9o) in the {x,y,z) system. 



where {xs,ys) are the angles of the source in the ob- 
server's sky. Inverting Eqs. (|23|l and (|24(l for the polar 
and azimuthal angles, it is convenient to consider good 
source alignments (e.g. small Xg and ys). The moti- 
vation for this approximation will become clear when we 
consider the magnifications of the relativistic images (de- 
fined below). If the observer is not exactly at = 0, we 
can approximate: 9s ~ tt — 9o — S9 and 4>s ~ tt — (50 with 
(56' <C 1 and Scf) 1 (an observer exactly at 6'o = will 
be considered in Sec. V). Expanding Eqs. H23|l and H24|) 
to second order in the perturbations, one finds that 
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x'^ sin 29o 



sinQp 



Vs COS Uo 



Us COS Uo 



(25) 
(26) 



Now we want to know how the restrictions in the pa- 
rameter space (A, r]) and in the lens geometry are re- 
flected in the possible values of the image coordinates 
(ai,/9i). We begin with the photon region discussed in 
Sec. II. This constraint correspond to a closed region 
in the observer's sky. It is useful to write its bound- 



ary as a parametric curve £,{^p)^ where ^ 



+ y and 



tan ip = y/x. Inserting Eqs. (O, (HHJ, (EH) and in 
the definition of (using the small angle approximations 
a ~ TqX, P ~ Toy), we can solve for r{(p) and obtain the 
desired parametric curve. When doing so, one encoun- 
ters a sixth order polynomial in r. The largest positive 
root is valid for —tt/2 < (fi < Tr/2 and the second largest 
positive root for tt/2 < tp < 3tt/2. 

Another constraint comes from the polar movement of 
the light rays. In Sec. II we pointed out that the motion 
in was restricted between the turning points 6'niin/inax 
defined by Eq. H15|l . Therefore, points in the parameters 
space where the inequality 0,nin < , < ^max is not 
satisfied must be discarded. 

Now we will prove that if A and rj are given by Eqs. 1)21(1 
and H22I) . then the inequality is always satisfied for 9o- 
Using the notation: uj = cos^Oj, Wj = sin^ 0j (j is any 
subscript), the above inequality is equivalent to u,n > 
(m stands for "min/max"). In this notation, the condi- 
tion for a turning point in 9, Q{9) = 0, becomes 



r]+ (a^ 



V 



2 2 

a M„ 



Writing 
get 



. (27) 
X and substituting into Eq. H27|l . we 



1 

2^ 



± 



[a^wo - («f + /3f)f - ia^lSfwo 



(28) 



On the other hand, using Eqs. (|21|l and H22f) we have 

+ (3f = a^uo + \^ + ri. (29) 

The "radius" A^ -I- 77 must be greater or equal than the 
boundary of the photon region defined by Eqs. I|17() and 
(|18|l and the line 77 = 0. It is then easy to show that the 
minimum value of this "radius" is reached when 77 = 
and a — 1. The actual minimum value is thus (A^(r) -|- 
^(r))min — 4. Therefore, by Eq. we have that of + 
Pi > 4, and since < a^Wo < 1, it follows the inequality 



a'^Wo — {af + Pf) < and hence, x < 0. The equality 
is satisfied only for /? = 0. We then can conclude that 

If Us > Uo, we can ensure that Um > Ug by doing the 
same construction in the frame of the source; 



\ 2 2 

A ^l^Ws , 



(30) 
(31) 



where ji and Qi are the coordinates in the source's sky. 
By using Eqs. ^ a.nd it follows that 



+ a'^{us-Uo), (32) 
(33) 



7^ 



a. 



Since for physically relevant situations we must have 
f2| > (where the equality holds only for Um — Us), 
the excluded points will be those between the curves 



A = ± 



Uo) 



1/2 



(34) 



where Us > Uo and of > a^Wg. The forbidden regions 
in the observer's sky are shown in Fig. Ofor the case of 
Us > Uo- For Us < Uo, only the photon region is present. 
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FIG. 3: Forbidden regions in the observer's sky for Us > Uo- 
The forbidden regions are shown in shaded gray. In this case 
we used a = 1, Uo = 1/2 and its = l.luo. For Us < Uo only 
the photon region is present. 

Next, we focus in relating the null equations of motion 
considered in Sec. II with the geometry of the gravita- 
tional lens system. In the familiar case of a Schwarzschild 
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black hole lens, the movement of light rays is restricted 
to a plane. Therefore, without loss of generality we can 
work at the equatorial plane: = 7r/2. Then, by substi- 
tuting Eq. 11U|) in and setting a = 0, we obtain the 
familiar Schwarzschild deflection angle: 

A</, = 2 r , (35) 

where 77 = by the fact that we are working at the 
equatorial plane [see Eq. Hlt)|) ]. This deflection angle can 
be written in terms Xi using Eq. (|21|l . Then, by use of 
the familiar "lens equation" we can solve for the position 
of the "virtual" images 0. However, another approach 
is to give A<^ in terms of the geometry of the observer- 
source pair and solve Eq. for Xi. In this case we 
have 

^c^=[-^r^v^ ^ '"^""n (36) 
^ I 27r(n + 1) - 0s , < ^ ' 

where n — 0, 1, 2 . . . is the number of windings around 
the z axis. The solutions for n = are the familiar 
weak field images, and for n > 1 we have the relativistic 
images studied in 0, |^ ^ . This last method is more 
useful when we consider trajectories outside the equa- 
torial plane. In this case, Eqs. and with H36() 
become our "lens equations" . Nevertheless, for trajecto- 
ries outside the equatorial plane, an additional parameter 
appears: the number of turning points in the polar coor- 
dinate (m) . For the case of the Schwarzschild black hole 
it is easy to prove that m is related to n as m = 2n -I- 1 
by the fact that photons travel only in a plane. 

What about the case of a Kerr black hole? In this situ- 
ation, since the movement of light rays is not necessarily 
restricted to a plane, we must consider n and m as inde- 
pendent parameters. Additionally, in trying to write A0 
as in Eq. one is faced with a possible complication: 
unlike in the Schwarzschild space-time, the Kerr geom- 
etry admits turning points in cj) which could complicate 
the analysis. By using Eq. ||2Jl with = is easy to 
show that the possible turning points occur at 

r±{0) = ^{(A-asin^^) 

±^{\-a sin' Of - a?X^ cos^ . (37) 

Both solutions are surfaces of revolution around the ro- 
tation axis of the black hole. It is easy to prove that for a 
given A, only one of the surfaces lies outside the horizon. 
Since in our case we are considering photons that come 
from infinity, reach a single turning point in r and return 
to infinity, they can cross such a surface at most in two 
points. In consequence, there can be at most two turning 
points in (j). However, by using Eqs. Q and 10} for an 
observer and a source at large r, it is easy to prove that: 
r1 sin' 9s{d(j)/dt)s « A « r' sin' 9oid(j)/dt)o (conservation 



of angular momentum). Therefore, the asymptotic sign 
of d(j)/dt must be unchanged by the gravitational inter- 
action. This rule out a single turning point in (j) which 
would obviously change that sign. We conclude that the 
number of turning point in <j) must be zero or two. In that 
case, the most general expression for is still given by 
Eq. (|35|) but with the following modification: the first 
case is to be use when the right hand side of Eq. is 
negative, and the second case otherwise. 

By using Eqs. (EIJ, (US, (123 and (HHIl, the "lens equa- 
tions" ifTUI) and lfTT|l can be expressed as 

H{xi,yi) - G{xi,yi,Xs,ys,m) = , (38) 

A<p{xs,ys,ri) - L{xi,yi) - J{xi,yi,Xs,ys,m) = , (39) 

where H and L [G and J) are the radial (angular) inte- 
grals of Eqs. (ITU)) and lfTT)l respectively, and n = 0, 1, 2, . . . 
and m — 0, 1, 2, . . . are the number of windings around 
the z axis and the number of turning points in the polar 
coordinate respectively. 

Although the integers n and m should be considered 
independent, as we will see in Sec. IV, the magnifica- 
tion does not depend directly on n. Therefore, in this 
paper we consider m to be the fundamental parameter. 
This interpretation is confirmed by our numerical results 
where we find that the magnification always decreases as 
we increase m, and for a given m, it can ever increase for 
images with larger n (see Sec. V). 

In the familiar Schwarzschild gravitational lens, we al- 
ways have two images which are formed by light rays 
which suffer small deviations in their trajectories. They 
are called weak field images. Additionally, we have an 
infinite set of faint images at both sides of the black 
hole jJl. These are relativistic images^ and are the re- 
sult of photons that orbit the black hole several times 
before reaching the observer. In Kerr space-time, where 
photons' trajectories are no longer confined to a plane, 
the concept of "orbiting the black hole several times" 
can be subtle. For this reason, we classify images as fol- 
lows: images with m — Q,l are called direct images (DI 
hereafter), and images with m > 2 are called relativistic 
images of order m (RI hereafter). 

As we show in the Appendix, Eqs. (|38|l and (|39|l can be 
written in terms of elliptic integrals, and are highly non- 
linear in all arguments with the exception of n and m. 
Therefore, their solution require numerical and graphical 
methods. To solve the lens equations, we use the stan- 
dard routine "FindRoot" build in MATHEMATICA 0. 
To give an initial approximation and to ensure that a so- 
lution exist at all, we use a variety of graphical methods. 
Since, as we will see in Sec. V, RIs appear very near the 
boundary of the photon region, we found that is useful 
to write: Xi — {£,(}p) -\- 5) cosip and j/i = {£,{^p) -\- 5) sini^, 
where < 1^9 < 27r, 5 > and ^((/s) is the boundary 
of the photon region as discussed earlier in this sec- 
tion. We then express all functions in terms of 5 and 
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Lp to avoid entering the photon region. We find that 
RIs form at very small 5. To find the images we plot 
the surfaces zi{5,Lp) = H{S,ip) — G{6,(p,Xs,ys,m) and 
Z2{5,(p) = A(f>{xs,ys,n) - L{S,(p) - J{S,ip,Xs,ys,m) for 
a given source position and image numbers (n,m). The 
intersection of both surfaces with the plane z = form 
various curves in that plane. If the curves formed by the 
two surfaces intersect each other, there is a solution. Vi- 
sual inspection of these curves allow us to give an initial 
approximation for the numerical routine. As an exam- 
ple, we plot part of both surfaces in Fig. 0] for n = 3 
and m — 4. We usually do this for the two lowest val- 



0.03 




0.02 



0.01 



FIG. 4: Here we plot the surfaces zi{S,ip) = H(5,ip) — 
G{5,ip,xs,ys,m) and Z2{5,ip) = A(t){xs,ys,n) — L{5,ip) — 
J{5,tp,Xs,ys,'rn) for a source at Xs — ys = l/\/2 arcmin, an 
observer a,t 9o = n /A and for n — 3, m = 4, a — 0.9999. We 
have chosen ro ~ rs — 8.5 kpc as in the numerical examples 
of section V. The intersection of the surfaces with the z — 
plane (gray) form various curves in that plane. The intersec- 
tion of two such curves marks the position of a solution to 
the "lens equations". The range of z has been restricted to 
—0.01 < z < 0.01 in order to clarify the intersection of the 
curves. 

ues of TO where a solution can be found. The possible 
values of n are bounded from above by the requirement 
\A(j>\ < maxlL^S^cp) + J{5,(p,Xs,ys,'m)\- Inconsequence, 
for a given m we usually find images for just a few values 
of n. This procedure is very tedious and time consuming, 
and for these reasons we are unable to give a complete 
phenomenological description of the behavior of the im- 
ages for any given observer-source geometry. Instead, we 
shall present some numerical examples of the kind of be- 
havior that can be expected by having a Kerr black hole 
as a gravitational deflector and null geodesic motion off 



the equatorial plane. 

As pointed out before, several authors have already 
obtained analytic approximations to the radial integrals 
in the "strong field limit" for Schwarzschild, Reissner- 
Norsdstr0n and Kerr black holes ^, J^, 18} . Their 
approximations are valid for relativistic photons, with 
Tmin close to T and for small deviations from the equato- 
rial plane in the case of the Kerr black hole. However, 
since in this article we are interested in orbits that can 
deviate significantly from the equatorial plane we can no 
longer use such approximation schemes. Moreover, their 
whole approach fail for the angular integrals because, in 
general, there is no relation between rmin and the turning 
points in 9. In Sec. V, to make contact with the relevant 
literature, we consider trajectories close to the equatorial 
plane to be compared with Bozza (18| procedures. 



IV. THE MAGNIFICATION IN A KERR 
GRAVITATIONAL LENS 

The magnification of an image is defined as the ratio of 
the observed flux to the flux of the unlensed source. By 
Liouville's theorem, the surface brightness is unchanged 
by the gravitational light deflection Therefore, the 
magnification is defined as the ratio of the solid angle 
subtended by the image to the solid angle of the unlensed 
source: 



duji 



1 

Jj\ 



(40) 



where J is the Jacobian of the transformation (xs,ys) 
{xi.yi). Writing Xs = Xs{xi,yi) and = ys{xi,yi) 
we can find expressions for dxs/dxi, dxs/dyi, dys/dxi 
and dys/dyi by differentiating the lens equations H38|) 
and (|39|l with respect to Xi and yi. After some algebraic 
manipulations, we find that 



M ■ 



(41) 



where 





dG 




dG 


ai ^ 


dxs 




dys 


d(j)s dJ 






dJ 


dxs dxs 


ai — 


dVs 


dys 


_ dH dG 




dH 


dG 


dxi dxi 


dVi 


dyi 


dL dJ 


, Pi 


dL 


dJ 


' dxi dxi 


dyi 


dyi 



(42) 

"3 = -TT— - TT— , "4 = -TT- - TT- 7 (43) 

(44) 

dyi ^ dyi ' ^^^'^ 

The above derivatives are very cumbersome and we will 
not expand it, neither show explicitly them here. In fact. 
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instead we shall use the equivalent numerical derivatives 
in all forthcoming calculations. Note that the magnifica- 
tion do not depend on n directly, since we have to take 
the derivative of A(/). Because of the complexity of these 
expressions, we are unable to give a complete description 
of the caustic structure of the Kerr space-time. Thus, we 
shall limit ourselves to compute the magnification of the 
images that we found. For a description of the caustic 
structure of a Kerr space-time, see |^ . 



V. RELATIVISTIC IMAGES BY A KERR 
BLACK HOLE GRAVITATIONAL LENS. 

In this section we present numerical calculations for 
the positions and magnifications of the RIs for different 
source-observer geometries. Our purpose is to provide a 
physical insight of the phenomenology that can be ex- 
pected from a rotating black hole behaving as a gravi- 
tational deflector. Also, our general procedures and nu- 
merical solutions can provide a set of "test-bed" calcula- 
tions for more sophisticated gravitational lens models to 
be developed in the near future. To be able to compare 
our results with recent published articles |^ 0, 0, 0| 
we shall consider a gravitational lens composed of a ro- 
tating black hole at the Galactic center with a mass of 
M = 2.8 X 10^ Mq, and located at a distance of Vo = 8.5 
kpc We shall also take Vg — Vo- 

The section is subdivided as follows: we first consider 
two simple cases involving a Schwarzschild black hole 
(a — 0) and an observer located at the pole {9o = 0). 
Next, as a consistency check, we consider an observer 
at the equator (9o — 7r/2 ) and compare our calculations 
with those of Bozza [ig- Finally, we work out a more 
general case of an observer at say 9o = tt/4. The RIs are 
classified as follows: we use a plus {+) sign for images 
with sign = sign ys (e.g. in the same "side" of the 
source) and a negative (— ) sign otherwise. We do all cal- 
culations for the two lowest values of m where a solution 
exist. We find that, in general, images with larger m are 
more demagnified. The following notation will also be 
used: ^ = 



' x'^ + , ip = arctan(2//a;), where {x,y) are 
the angular positions in the observer's sky. 

In all geometric configurations considered in this sec- 
tion, we found that for to = 1 we recovered the usual 
weak field images without any noticeable effect from the 
black hole's spin. Moreover, we could not find any image 
with 771 = for this kind of geometry. 



A. A Schwarzschild Black Hole as a Gravitational 
Deflector 

In the Appendix we show that for a — 0, the angu- 
lar integrals can be solved in close form. As expected, 
we find that the RIs are always found in the line join- 
ing the source position (xg^ys) and origin of the ob- 
server's sky regardless of the inclination of the observer 



9o- The separation of the two outermost images is 
always the same regardless of the source position (or 
9o)- They are: ~ 16.952 /^arcsec (777 — 3) and 

^_|. = ^_ « 16.931 ^arcsec (777 = 5) for the two lowest or- 
der images. These results are close to those found in P]: 
^+ = ^_ w 16.898 /zarcsec and = ^_ w 16.877 ^arcsec 
respectively. The magnifications calculated by Eq. H41|) 
are also close to those in the literature. However, we en- 
countered problems of numerical noise when using the 
exact expression for J(xi, Zg, j/s, to). For that rea- 
son, to check that we obtain the right magnifications, 
we used the expression for a ^ [Eq. (|A.38|l ] and set 
a ~ 10""'^^. The magnifications obtained this way for the 
outermost image as a function of the source separation 
agree with the literature Q. For instance, we obtained 
/i+ = M- « 3.5 X 10-18 and Ai+ = « 3.6 x IQ-^^ 
for = 1 arcsec and £,s = 100 /iarcsec respectively. For 
comparision, in Ref. (sj we read: fi+ = fj,- ^ 3.5 x 10^^^ 
and fi+ = n- !v 3.5 x 10~^'^ respectively. 



B. An Observer at the Kerr Black Hole's Pole 

(do = 0) 

Although is very unlikely that an observer will be ex- 
actly at 9o = 0, this is the simplest case that can be 
solved quasi-analytically. Using Eq. H21|) for 6*0 = we 
find that A w 0. The minimum value of rj, corresponding 
to the photon sphere, is found from Eqs. H17|l and H18() 
by setting A = 0. We obtain 



'7min(a) 



r(a)^ — 



(46) 



where 



7:(a) = -(3-a2)[i + 2cos(VV3)] , (47) 



tani/j 



3^/3aVl08 - 135a2 + 36a4 - 4a6 
54 - 81a2 + 18a4 - 2a^ 



(48) 



Eq. H46|l is a decreasing function of a satisfying 14.33 < 
Vmin ^ 27. Since we can relate the position of the images 
to r] by Eq. d « {l/ro)^/v + a^, Eq. lO sets a 

lower bound to the separation of the images d. 

The null equations of motion simplify considerable for 
A = 0. First we note that from Eq. (6), there are no 
turning points in 9 since 77 > 0. Therefore, the angular 
integration limits can be shown to be 



P 



+ 2(277 + 1) 



7r/2 



(49) 



where p = ±1 is the parity of the image ( — for images 
that are in the same side of the source when a = and + 
for images in the opposite side). Here, 77 is the number 
of loops around the black hole. To relate 9s to the un- 
perturbed position of the source in the observer's sky we 
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can no longer use the approximate expressions H25|) and 
(I26|l . Instead, it can be shown that 



(50) 



- 6 1 



where we are assuming that + fo/Ts) ^ 1- Now, 
using Eqs. (|A.1|I . 1A.4|I and l|A.30() given in the Appendix, 
the null equations of motion (fTH)) and ltTT|l become: 



F{i;,k) 



\/{ra - rc){rb ~ r^) 
1 



pF 



6 1 



-2{2n+l)K 



= 0, (51) 



Aip = 



4a 



V(ra-rc){rb-rd) ^ 



-/3,F(^,fc)] , 



where 



arcsm 



rd 



in - rc){ra - rd) 
{Ta - rc){rb - Td) 

jrj - n){ra - rd) 
(rj - ra){rb - rd) 



VT^+ (-1)^+1 



VT~^ ((-l)Vl-a2 +ra-l 



r, = l + (-l)J+Vl-a" 



(52) 

(53) 

(54) 

(55) 

(56) 

(57) 
(58) 



and r„ = 



''b, ?'c, r.i {ra > n > rc > ra) are the roots 



of R{r) = 0, and Eq. H51|) is valid only for a 1 (for a = 1 
the right hand side of Eq. (|5^ has to be substituted by 
Eq. HA.11() of the Appendix) . The procedure to calculate 
the positions of the RIs is the following: for a given source 
separation ^s, we use Eq. H51|) to calculate the angular 
separation of the RI (^i) and then insert this value in 
Eq. 1)52(1 to obtain the offset from the source inclination 
Aif = (pi — (ps. In calculating one has to use Eq. (|46|l to 
determine its minimum value and avoid numerical prob- 
lems. Note that for a = 0, Eq. (|FT|l gives ipi = ips as 
expected. In Fig. O we plot A(p for the three outermost 
images (n = 1, 2, 3) as a function of a for a source located 
at = 1 /iarcsec. The sign of A(p is that of a (which is 
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FIG. 5: Deflection in the image inclination {Aip = ips) as 
a function of the normalized black hole angular momentum 
a. Here we consider an observer located at = and a 
source located at = 1 /iarcsec. The parameter n is the 
number of loops around the black hole. These plots are for 
p — —1, but photons with p = 1 have almost the same curves 
(indistinguishable within the plot resolution). 



what one intuitively expects). For a given n, both images 
with opposite parity have almost the same Aip. However, 
the greater the number of loops n, the greater the deflec- 
tion Aip for a given a. The physical picture emerging is 
a very simple one: the angular momentum of the black 
hole just adds a "twist" in the direction of the rotation 
to the usual Schwarzschild trajectory. To give an intu- 
ition of the full movement of the images, including their 
separation , in Fig. |H1 we have plotted the position of 
the two lowest order images (rt = 1,2) as a function of 
the spin parameter a. 

What about the magnifications? In this simple case 
where we have circular symmetry as seen by the observer 
the magnification is given by 
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d9_ 



dg 



(59) 



where g — g(Cj, Cs) is the left hand side of Eq. H51|l . This 
expression for the magnification can be verified by cal- 
culating the Jacobian of the transformation (6,<i5s) 
(Cii 'Pi)- We find that, as in the Schwarzschild case, both 
images with the same winding number n have approxi- 
mately the same magnification (to at least three signif- 
icant figures). In Fig. [3 we plot the magnifications for 
the three lowest order images (n — 1,2,3) and with a 
source located at = 1 //arcsec. The net effect of the 
angular momentum in this case is to enhance the bright- 
ness of the images. For the Schwarzschild case (a = 0), 
the magnifications obtained with Eq. H59|) agree perfectly 
with the results of Ref. ^ . 
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FIG. 6: Positions of the two lowest order primary RIs (n — 
1,2 and p = —1) as a function of the normalized angular 
momentum of the black hole a. The arrows in the curves 
represent the direction of the movement as we increase a from 
to 1. The source is located at = 1 /larcsec and ips = 7r/4. 
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FIG. 7: Magnification for the three lowest order RIs (n — 
1, 2, 3) for an observer at = and a source at = 1/^arcsec. 
Here we consider the case of p = — 1, but the magnifications 
for p — 1 are almost the same (indistinguishable in this plot). 



C. An Observer at Kerr Black Hole's Equator 

(0o = 7r/2) 

In this subsection we consider an observer located ex- 
actly at 9o = 7r/2. Our purpose is not to give a com- 
plete account of the phenomenology of the RIs since, 
from an astronomical perspective, it is very unlikely that 
an observer will be exactly at the equator (also setting 



6o = 7r/2 does not significantly simplify our analysis). 
Rather than that, this special case will serve as a con- 
sistency check by allowing us to compare some of our 
numerical results with those of Bozza [iSj . 

In that article, the author considered "quasi- 
equatorial" orbits in Kerr spacetime. One of the main ap- 
proximations employed was that the horizontal position 
of the RIs {xi in our notation) was calculated indepen- 
dently using the familiar lens equation in the equatorial 
plane. This, of course, assumed that the motion in (jj was 
unaffected by the motion in 9. Therefore, we are forced 
to consider very small source declinations (ps ^ I. To 
this end (and for simplicity) we fix the source at = 1 
arcsec and ys = I /larcsec {(ps ~ 10^^). Then, we cal- 
culate the position and magnification of the lowest order 
RI in the same side of the source as a function of the 
spin parameter a. We fixed n = 1 to be consistent with 
the reference, where the RIs are classified only with the 
winding number n. We also use a < for latter con- 
venience. The results using the approximate equations 
(66), (77) and (90) of Ref. [lff| and the equations of this 
paper are shown in Table I [26]. 

The general behavior of the position of this image for 
a < is the following: as we increase |a|, it moves toward 
the X axis and gets very close to the forbidden region 
defined by Eq. Then for \a\ > 0.7 it starts to 

increase ipi until a maximum value of ipi « O.Itt is reached 
(when \a\ —>■ 1). This increase in ipi coincides with a 
change in the order of the image around a ^ 0.6 — 0.7: 
from m = 3 to m = 2. This illustrates the importance 
of considering m as an independent parameter. Also, it 
gives a better physical intuition of the qualitative form 
of the photon's trajectory (see Fig. ISJ. 



source 




obssrver 



FIG. 8: A typical light ray trajectory with m = 2 and n = 1 
for an observer at the equator [6o ~ 7r/2). The cones rep- 



resent the turning point angles 6, 



mill / max • 



The number of 



turning points in a trajectory is the number of times the pho- 
ton touches these cones (in this case, two). The drawing is 
not to scale. 

The calculations done with the approximate expres- 
sions of Bozza [l^l reproduce the behavior of this RI 
very well until \a\ gets higher than \a\ ~ 0.9 and the image 
moves too far away from the equator where the "quasi- 
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TABLE I: Here we fix the source at = 1 arcsec, = 7r/4 and calculate the position and magnification of the lowest order 
RI in the same quadrant as the source as a function of the black hole normalized angular momentum a. The numerical results 
using the approximate expressions of Ref. fl^ are shown in columns 2-4. Our results are shown in columns 5-8. We also 
show the number of turning points for each case (m). All angular positions are given in microarcseconds and all angles in 
radians. 
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a 


6 




TV 






M 
















m 


-10"" 


16.952 


2.615 X 10 


— 


2.9 


X 


10" 




16.952 


2.618 X 10' 


— 


2.9 


X 


10" 


"TH~ 


3 


-0.1 


16.295 


1.442 X 




11 


1.8 


X 


10" 


22 


16.296 


1.444 X 10" 


11 


1.8 


X 


10" 


22 


3 


-0.2 


15.619 


7.161 X 


10" 


12 


1.0 


X 


10" 


22 


15.619 


7.174 X 10" 


12 


1.1 


X 


10" 


22 


3 


-0.3 


14.917 


4.832 X 


10" 


12 


8.1 


X 


10" 


23 


14.918 


4.844 X 10" 


12 


8.5 


X 


10" 


23 


3 


-0.4 


14.187 


3.758 X 


10" 


12 


7.4 


X 


10" 


23 


14.188 


3.771 X 10" 


12 


7.6 


X 


10" 


23 


3 


-0.5 


13.422 


3.223 X 


10" 


12 


7.5 


X 


10" 


23 


13.423 


3.240 X 10" 


12 


8.1 


X 


10" 


23 


3 


-0.6 


12.612 


3.024 X 


10" 


12 


8.5 


X 


10" 


23 


12.614 


3.050 X 10" 


12 


1.8 


X 


10" 


22 


3 


-0.7 


11.743 


3.155 X 


10" 


12 


1.1 


X 


10" 


22 


11.747 


3.207 X 10" 


12 


1.2 


X 


10" 


22 


2 


-0.8 


10.790 


3.881 X 


10" 


12 


1.6 


X 


10" 


22 


10.802 


4.051 X 10" 


12 


1.7 


X 


10" 


22 


2 


-0.9 


9.695 


6.672 X 


10" 


12 


3.3 


X 


10" 


22 


9.741 


8.428 X 10" 


12 


5.0 


X 


10" 


22 


2 


-0.95 


9.018 


9.763 X 


10" 


12 


4.6 


X 


10" 


22 


9.147 


3.672 X 10" 


11 


2.8 


X 


10" 


21 


2 


-0.99 


8.063 


4.058 X 


10" 


12 


9.7 


X 


10" 


23 


9.042 


0.1167 




7.0 


X 


10" 


22 


2 



using the formulas of this paper 



equatorial" approximation is no longer valid. Choosing 
a > for a source in that position would increase ipi 
from the beginning and would not allow us to do the 
comparison with the reference. 

For the magnifications, we found that our results 
roughly agree with those of Bozza, except for a > 0.9. 
However, we have to admit that we encountered prob- 
lems with the numerical derivatives for large values of a. 
Therefore, our results for this case should serve only as a 
order-of-magnitude estimate. Nevertheless, both results 
agree in that the net effect of the black hole spin is to fur- 
ther demagnify the RI as compared to the Schwarzschild 
case (a = 0). 



D. An Observer Off the Kerr Black Hole's 
Equatorial Plane {9o = 7r/4) 

This is a more general case that will allow us to illus- 
trate the phenomenology of the RIs in a more realistic 
astronomical situation. We found that the general be- 
havior of the RIs for 6'o — 7r/4 is valid for a wide range of 
inclinations (but the precise location and magnifications 
of the images will depend on 9o)- 

We start with a slowly rotating black hole and set 
a = 10~^. We then calculate the positions and magnifi- 
cations of the two lowest order RIs as a function of the 
source separation. More precisely, we consider a source 
travelling along the dashed line shown in Fig. O In the 
same figure, we show the positions of the two brightest 
images which correspond to m = 3, n = 1 in all cases. 
The numerical results for the two lowest order RI found 
are shown in Table II for m — S and Table III for m — b. 

Both images appear in opposite sides of the black hole 
as in the Schwarzschild case. However, for the slowly 
spinning black hole, they rotate counterclockwise as the 
source approaches the origin of the observer's sky and 
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FIG. 9: Position of the lowest order RIs (m = 3) as a function 
of the source separation. The black hole spin parameter has 
been fixed to a = 10~". The source is travelling along the 
dashed line starting at Xs = i/s ~ 1 arcsec and finishing at 
Xs = Us — —1 arcsec. The arrows illustrate the direction of 
movement of the two RIs, and their starting positions have 
been label by the word "start" . The solid circle is the bound- 
ary of the photon region. 

emerges in the opposite quadrant (see Fig.l^)). If we take 
a < 0, the effect is the same but the images rotate clock- 
wise. Note that for very large source separations, the 
positions and magnifications of the images tend to the 
familiar Schwarzschild values |8| . 
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TABLE II: Positions and magnifications of the lowest order RIs (m = 3) as a function of the source separation for a slowly 
rotating black hole. The spin parameter of the black hole is fixed at a = 10~® and the source is travelling along the dashed 
line of Fig.|5| Here a positive value of ^3 means that the source is in the first quadrant of the observer's sky {ips — 7r/4) and a 
negative value means that it is located in the third quadrant {ifs = 37r/4). In this case we found solutions only for n = 1. All 
angular positions are given in microarcseconds. 
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TABLE III: Same as Table II but with m — 5. In this case we found solutions only for n = 2. 
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2.0 X 


10" 


20 


-10^ 


-16.931 


-3.6657 X 10"^ 


2.0 X 


10" 


20 


16.931 


3.6657 X 10"^ 


2.0 X 


10" 


20 


-10* 


-16.927 


-0.35624 


1.9 X 


10" 


20 


16.927 


0.35624 


1.9 X 


10" 


20 


-lO'^ 


-16.672 


-2.9498 


1.6 X 


10" 


20 


16.672 


2.9498 


1.6 X 


10" 


20 


-lO'' 


-13.984 


-9.5443 


5.2 X 


10" 


21 


13.984 


9.5443 


5.2 X 


10" 


21 



If we compare the magnifications with those calculated 
in Ref. @ for a source at the same separations but with 
a = (Schwarzschild black hole), we again observe that 
the effect of the black hole spin is to further demagnify 
the RIs. For instance, for a source at = 1/J,arcsec, the 
magnifications calculated in 8] for the two brightest RIs 
were: /i+ = /i_ w 3.5 x 10~^^ while our calculations for 
a = 10~^ yield tJL+ = H- k. 2.3 x 10"-^'^, a difference of 
five orders of magnitude! A similar behavior is observed 
for m = 5 (see Table III). This example illustrates how 
even a small value of a can make significant differences in 
the phenomenology of the strong-field gravitational lens. 

Now we consider a fixed source at = 1 arcsec, 
Lfs — -K /A and calculate the positions and magnifications 
of the two lowest order RIs as a function of the spin pa- 
rameter a. We consider the case of a > only for simplic- 
ity. To illustrate the phenomenology of the RIs we have 
included Fig.ITOI We observe that forl0-*^<a<10-3 
the two lowest order sets of RIs have m = 3 and m — 5 
respectively as in the case of a Schwarzschild black hole. 
We also observe two images per value of to. As we in- 
crease a, these images move toward the x axis of the ob- 
server's sky. For a > 10~^ these images disappear in the 
X axis and give rise to images of order to = 2 and to = 4. 
For TO = 2, we observe only one image in the third quad- 
rant of the observer's sky. On the other hand, for m = 4 



we found two images (one in the first quadrant and one 
in the third). Increasing a even more (a = 0.9999) we 
discovered that for to = 2, two images are formed in the 
third quadrant while for m = 4, a total of six images 
are observed: two in the first quadrant and four in the 
third quadrant. The positions of the lowest order images 
(to — 2 and to = 3) are shown in Fig. 1101 while the posi- 
tion and magnifications of all images (to = 2,3,4,5) are 
listed in Tables IV - VII. Is interesting how the number 
and order of the images depend on the value of the spin 
parameter a. 

With regard to the magnifications, we observe that as 
in previous cases, the effect of the black hole angular mo- 
mentum is to demagnify the images (although there is a 
slight increase in the magnification for a = 0.9999). Also 
note that the images with to = 4 in the first quadrant 
have magnifications very similar to those with to = 2 
(located in the third quadrant). Therefore, these im- 
ages should be considered "dual" to each other since they 
form the pair of brightest images for a given value of a 
(a > 10^^). An interesting consequence of the black hole 
rotation is that for large values of a (a > 0.1) the two 
brightest images have different magnifications. This is to 
be compared to the case of a Schwarzschild black hole 
where the two brightest RIs have exactly the same mag- 
nification. Also, for large a the RIs are very static as 
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TABLE IV: Positions and magnification of the order m = 3 Rls as a function of the black hole spin parameter o. Here we 
consider a fixed source at = 1 arcsec and ips = 7r/4. For a > 10"^ these images disappear in the x axis of the observer's sky. 
All images have n = 1. All angular positions are given in microarcseconds. 



a 






/' 








/' 




10-" 


10.461 


13.340 


3.9 X 10" 


iK 


-10.4605 


-14.930 


3.9 X 10" 


18 


lo--' 


-12.832 


11.077 


3.2 X 10" 


18 


12.8318 


-11.0777 


3.2 X 10" 


18 


10-* 


-16.932 


0.82137 


2.4 X 10" 


19 


16.9326 


-0.82145 


2.4 X 10" 


19 


10"^ 


-16.947 0.066386 


2.3 X 10" 


20 


16.9565 


-0.0664584 


2.3 X 10" 


20 



TABLE V: Same as Table IV but with m = 5. All images have n = 2. 



a 


4 


yt 








Vi 






IQ-^' 


7.9844 


14.930 


8.1 X 10" 


21 


-7.9844 


-14.930 


8.1 X 10" 


21 


IQ-^ 


-16.330 


4.4712 


2.4 X 10" 


21 


16.330 


-4.4713 


2.4 X 10" 


21 


IQ-* 


-16.926 


0.36919 


2.0 X 10" 


22 


16.927 


-0.36923 


2.0 X 10" 


22 


IQ-^ 


-16.926 


9.5138 X 10"* 


2.0 X 10" 


23 


16.935 


-9.5469 X 10"^ 


2.0 X 10" 


23 



20 
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FIG. 10: In this plot we show the position of the lowest order 
Rls as a function of the normalized angular momentum a. The 
source is fixed at = 1 arcsec, ips = 7r/4. While we increase 
the angular momentum starting from a = 10~®, two images 
with m = 3 (boxes, ■) move toward the x axis as indicated 
by the solid arrows. Increasing a further (a > 10"^), these 
images disappear in the x axis and give rise to a single m = 2 
image (stars, -k) that move as indicated by the dashed arrow. 
For a = 0.9999 we find two images with m = 2 in the third 
quadrant of the observer's sky (crosses, x). 



in the Schwarzschild case (although for a = they have 
'Pi = 'Ps)- In other words, their positions do not depend 
on the location of the source. 

What is even more surprising is that for a ^ and 
within our numerical precision, it seems that the ratio of 



the magnifications of the two brightest Rls is insensitive 
to the position of the source. This raises the possibil- 
ity of extracting information about the orientation and 
spin of the black hole by comparing the brightness of this 
two images. However, to move this proposal forward we 
would need to study the behavior of the ratio of their 
magnifications for different values of 9o- We do not at- 
tempt to carry such analysis this time. We hope that any 
new approximation scheme developed in the future can 
allow us to address that question. 



VI. CONCLUSIONS 

In this article we have explored the phenomenology of 
strong field gravitational lensing by a Kerr black hole. In 
particular we have developed a general procedure to cal- 
culate the positions and magnification of all images for 
an observer and source far away from the black hole and 
at arbitrary inclinations. We have applied our developed 
procedure to the case of a black hole at the Galactic cen- 
ter with mass M = 2.8 x 10^ Mq and at a coordinate 
distance of Tq = 8.5 kpc. We have reproduced the posi- 
tions and magnifications of the lowest order relativistic 
images found in the references for a Schwarzschild black 
hole and for "quasi-equatorial" trajectories around a Kerr 
black hole. We have also presented new numerical results 
for the case of an observer located at 0o = 7r/4. 

Although we have not been able to give a full account 
of the phenomenology for all possible combination of the 
source-observer geometry and the spin parameter o, our 
limited results were useful to get a physical insight of 
the effects of the black hole angular momentum in the 
strong-field regime of gravitational lensing. Moreover, 
our findings can serve as a "test-bed" calculation for any 
improved new model of gravitational lens. 

There is no doubt that observations of the strong-field 
regime of gravitational lensing will be an extremely chal- 
lenging task in the near future. This is because, as we 
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TABLE VI: Positions and magnifications of the order m = 2 images as a function of tlie spin parameter a. Here we consider 
a fixed source at = 1 arcsec and (fs — tv/A. Tfiese images appear for a > 10~^. Note tliat for the first three values of a we 
have only one image but for a = 0.9999 there are two images. Also, in this case we found images with different values of n 
{n — 1,2). All angular positions are given in microarcseconds. 















a 




Vi 


M 




n 


10"^ 


-16.906 


-0.11746 


2.3 X 10" 


•21 


1 


0.1 


-16.439 


-1.2457 


2.5 X 10" 


22 


1 


0.5 


-13.200 


-6.1999 


8.3 X 10" 


23 


1 


0.9999 


-5.9311 


-11.907 


1.1 X 10" 


22 


1 




-9.6208 


-2.0417 


1.1 X 10" 


22 


2 



TABLE VIL Same as Table VI but with m = 4. Note that in this case we have two images for the first three values of a, but 
six images for a = 0.9999. Also, we find images with n = 1 — 7. 



a 


x+ J/+ fi+ n 


Vi fi n 


10'^ 
0.1 
0.5 


16.997 0.11757 2.3 x 10"^' 1 
17.354 1.2638 2.1 x 10"^^ 1 
17.680 6.6724 3.8 x 10"^^ 1 


-16.883 -0.26620 2.0 x 10"^" 2 
-16.244 -2.6972 2.6 x lO"^'^ 2 
-7.7361 -13.392 2.4 x lO"^'' 2 


0.9999 


9.1375 15.830 8.1 x lO"^'' 1 
3.4194 16.348 7.5 x 10"^^ 1 


-6.6670 -10.505 6.5 x lO"^'' 3 
-8.5479 -6.5516 8.2 x 10"^^ 4 
-9.0293 -4.2960 1.1 x lO"^"* 5 
-9.1975 -2.6233 1.7 x 10"^* 6 
-9.2611 -0.93129 6.5 x lO"^"* 7 



have seen, the relativistic images are always highly de- 
magnified. However, if we are able to observe them in any 
foreseeable future, they will provide one of the best tests 
of Einstein's General Theory of Relativity in strong grav- 
itational fields. Moreover, as we have seen, they could 
provide new tools to astrophysics by allowing the mea- 
surement of the orientation and/or magnitude of the an- 
gular momentum of the black hole. However, to fully 
confirm that this is the case, more research is needed to- 
ward developing an analytical solution of the strong-field 
gravitational lens problem in Kerr space-time. 



Acknowledgments 

S. E. Vazquez is very grateful to E. P. Esteban for all 
his support and advise during his undergraduate years. 
He would also like to thank NSF for a Graduate Research 
Fellowship and the University of California at Santa Bar- 
bara for a Broida Excellence Fellowship. E. P. Esteban 
thanks the support given by UPR-Humacao and Rice 
university during his sabbatical leave. 



APPENDIX 

Here we show how to write Eqs. Hl()|l and Hll|l in terms 
of elliptic integrals. All derivations are done following 
Ref. lii. 



1. Radial integrals 

We start with the radial integral of Eq. (|10|l : 



dr 



VW) 

4 



\/{ra - rc){rb - ra) 



F{i^,k), (A.l) 



where F{ip, k) is the normal elliptic integral of the first 
kind. Also, 



ip = arcsin 



n - Td 

Ta - 



^2 ^ {rg - rd){rb - Vg) 
{ra - rc){rb - ra) 



(A.2) 



(A.3) 



where Va = r,nin, n, r^, (r^ > n > > rj.) are the 
roots of R{r) = 0. For the radial integral of Eq. (|ll|l we 
have: 



L(A,?7) = 2 / — F=^dr 



1=1 "''^min 



dr 



{r-r,)^R(r) 
gY^T, [(l-/32)H(^,a2,fc)+/32F(^,fc) 



(A.4) 
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where !!(?/;, a^, k) is the normal eUiptic integral of the These expressions are only valid for a ^ 1. For a = 1 we 
third kind and have: 



-i + (-ir+ivr~ 



K, = l + (-1) 



»+il-«A/2 



a. 



ra - n 

in - n){ra ~ rd) 
in - ra){rb - rd) 
n - r„ 



n - n 



4a 



\/{ra - rc){n - Td) 



(A.5) 
(A.6) 
(A.7) 
(A.8) 
(A.9) 
(A.IO) 



2- A 2 
(r- 1)2 ^V~^ 

p. o , /3'(2-A) 



F(VAfc) + 2(1-/32) 



/32(2-A) 



r-a-l 



Ta-l 



(A.ll) 



where 



2(a2_i)(p_a2) 



a^E{u) + (fc2 - a2)u + (2a2fc2 + 2a 



3fc2)^(V',a^A:) 



a snu cnw dnu 



1 



(A.12) 



u = F{^,k), (A. 13) 2. Angular integrals 

2 (1 — rb)(ra — rd) . We begin by defining u = cos^9, such that 

a = r-, r , (A-14) 

(1 - rajin - Td) 1 (In, 

smed0 = — sign(7r/2-e')^ . (A.16) 

/^^~- (A.15) . . ' . ^ . 

-'^ Since all angular integrals considered in this paper are 

symmetric around 9 = 7r/2, we can restrict our inter- 
val to7r/2<0<7r and discard the sign of Eq. fulfill 
(remember that all integrals are positive definite). How- 
ever, we need to compensate for the case of a photon that 
crosses the equator. To this end we define the operator 

6*1 * 6I2 = cos 6ii cos 6*2 , (A. 17) 

Here E{u) is the complete elliptic integral of the second for any two angles 6*1 and ^2- Is clear that 9i * 02 > 
kind, and sn, cn, dn are the Jacobian elliptic functions. if both angles are at the same hemisphere and negative 
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otherwise. Therefore, is easy to show that we can write 

1*62 f~ fit 
± / = / + [1 - sign(0i * (?2)] / , (A.18) 

where u = min(Mi. U2) and u — max(Mi, 1*2). 

In the trajectories that we consider in this article the 
photon encounters various turning points defined by: 
Q{9) — 0. The positive solution to this equation is [see 
Eq. ©]: 



in 



1 - 



-I 1/2^ 



(A.19) 



where the subscript "m" stands for "min/max". The 
corresponding angles are 



^min / max — arCCOS (zb-^Mm) 



(A.20) 



where the plus sign correspond to ^min and the negative 

to ^niax- 

Now, for a trajectory that encounters m turning points 
(to > 1) we have 



± 



± 



± 



,i„/. 



...± 



m — 1 times 

+ [1 - sign(6'5 * 9,ns)] 
+ [1 - sign(6'o * e^no)] 
+2(to- 1) 



where 



arccos [sign{yi)^/u^] 
m odd : 



(A.21) 

(A.22) 
(A.23) 



with yi as the (possible) position of the image. In deriv- 
ing Eq. (TOTll we have use the fact that, as discussed in 
section III, u„i > Us, Uo- For to = we can write 



± 



[1 - sign(0, * Oo 



[1 - sign(^, * Oo)] 



where u = VDim(us,Uo) and u = max(us,Uo). 

To write the angular integrals of Eq. H10(l as elliptic 
integrals we change variables to u so that, for 9j and 
^min/max in the Same hemisphere, we have 



du 



±^/W) 2\a\ J^^ ^u{u,n~u){u-U3) 



hF{^j,K) 



where 



h = 



1 



arcsin Jl — 

Urn. 



Urn - U3 



(A.25) 

(A.26) 

(A.27) 
(A.28) 



U3 



1 f ^ + 7? 



1 - 



44 



1/2^ 



, (A.29) 



and Um > U3 (u3 < 0). On the other hand we have 



1 



du 



2|a| Jo \/u{u„i - u){u- Us) 
where 



hF{<^j,K) , (A.30) 



— arcsin < 



luj{Ur^ 



U3) 



Um{Uj ~ U3) 



(A.31) 



Using Eqs. lfOT)l . l|X24|) . ifOS)) and (fO0)l we can 

write the left hand side of Eq. (|10|l as 

G(A,77,0„to) = h{F{^s,^i) + F{^o,n) 

+ [1 - sign(0, * 0™,)] F($„ k) 
+ [1 - sign(0o * F{^o, k) 
+2(to- l)is:(K)} , (A.32) 

for TO > 1. Here, K{k) is the compete elliptic integral of 
the first kind: K{k) = F(7r/2, k). For to = we have 

Gi\,rj,e„m) = h{F{:^,K)- F{W,k) 

+ [1 - sign(0, * Oo)] F($, k)} , 

(A.33) 



where and ^ are given by Eq. (|A.27p with the substi- 
tutions Uj —^ min(us,Mo) an d Uj — > max(u<;,Uo) respec- 
tively, and $ is given by Eq. ljA.3ip with the substitution 



(A. 24) Mj — > min(ws. Mo). 
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Now we turn our attention to the angular integrals of where 
Eq. 1)11(1 . Making the usual change of variable, we get: 



,i„ / r 



de 



1 /■"'" du 
2|a| 7«j (1 - u)y/u{um -u){u~ M3) 

-n(v|/„7^«), 



1 — u. 



(A.34) 



where 



7 



Urn - 1 

The second integral we need is 

1 n du 



(A.35) 



2\a\ Jq (1 - U)yju{u„i - U){U - U3) 

i — U3 



Mm(l - M3) 
Urn - U3 



Using (jOU, C04)l . C04jl and 

(A. 36) integrals of Eq. (|lll) become 



(A.37) 



the angular 



J(A,??,0„m) = [n(*„72,K) + n(*„,72,K)] +-A^{(l-sign(0,*^?™,)) [F{<^,,k) ~ u^IV{<^s, .n) 

1 W-TD, -L 



"3 

+ (1 - sign(0o * dmo)) [-F($o, k) - U3n($o, k)] + 2(m - 1) [K{k) - ^/3^(p^ k)] } 

I 



(A.38) 



for m > 1, and 



J(A,77,0.,m) = [n(£,72,K)-n(vI/,72,«;)] +-^{[l-sign(0,*0„)] [F k) - u^n[^, p\ n)]] 

i — 1 — U3 



(A.39) 



for m = 0. 

For the case of a Schwarzschild black hole (a — 0) 
the angular integrals can be solved in close form since 
the polynomial in the square root become one of second 
order: 



where 



A2 + 77 



Similarly, 



(A.41) 



,in/r 



de 



du 



arcsin 1 — 2 



, (A.40) 



du 



2 Jo y/rju - {X'^ + r])i 







arcsin 


2 





(A.42) 



Therefore, using Eqs. I|A.21|) and (|A.24p . the right hand 
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side of Eq. (|10|l becomes 



for m > 1, and 
G{\r],6s,rn) = 



2^^ 



[27rm 



+sign(0s * 9„is) arcsin 1 — 2 — 



+sign(6'o * 6'mo) arcsin ( 1 — 2—2- 



(A.43) 



, arcsin 1 — 2- 
2V^ 1 2 V 



-siffn * I 



TT "Ml 

arcsm I 1 — 2 — 

2 V Um 



introduce a new variable: w = sin 9 = 1 — u. Thus, 



0, isin^ 61^6(61) 



1 

2 

1 

2W 



TT 



2A2u2 ^ 



„ , arcsin 1 — , , 
2 I - ^^,0 



(A.45) 



Similarly, for the — s- Uj integration we have 



dw 



2 Jn,, Wy/-X^ + (2A2 + T])W - (A2 + 77)m;2 



(A.44) 



1 

21^ 



arcsm 1 ; ^ 

2 I ^(1 - %) 



(A.46) 



for m = 0. 

The angular integrals of Eq. (|ll|l can also be calcu- 
lated in close form for a — 0. For that, is convenient to Therefore, the angular integrals of Eq. Hll(l become 

I 



J{\r],9s,m) 



sign(A) 



27rm -I- sign(0s * dms) arcsin 1 



77(1 - Us) 



for m > 1, and 
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arcsin 1 



2A2w 



77(1 - Uo) 



sign( 



) arcsin 1 



77(1 - Uo) 



(A.47) 



arcsin 1 



2A2?i- 



2, ,2 



77(1 - Mo) 



(A.48) 



for m = 0. 

Using the expressions for H , G, L and J given above, 
the "lens equations" (|10|) and pil) become: 

iI(A,77)-G(A,?7,0„77i) = 0, (A.49) 



where n = 0, 1, 2 . . . is the number of windings around 
the z axis, and A, 77 and 9s can be written in terms the 
observer's sky coordinates by using Eqs. l(?T|l. 
and lI23Il. 



A^-L{X,r])~J{X,r],9s,m) = 0, (A.50) 



with 



-(t>s - 2Tin , L + J < 
27r(72 + l)-(?!)5 ,L + J>0 



(A.51) 
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